Phase transitions in the antiferroelectric-ferroelectric mixed systems 
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Phase diagram of solid solutions of antiferroelectric and ferroelectric materials like PbZri-ajTi^^Os 
(PZT), Rb^(NH4)i-xH2P04 (RADP) is calculated at the region x < 0.1. Antiferroelectric and fer- 
roelectric order parameters (which correspond respectively to the difference and sum of two sublat- 
tices ions shifts) were introduced for antiferroelectric host lattice. Small admixture of ferroelectric 
component was considered as a random electric field source. The influence of this field on both 
aforementioned order parameters was calculated. Self consistent averaging with calculated random 
field distribution function made it possible to obtain the system of equations for order parameters 
of mixed system. The antiferroelectric Ta and ferroelectric Tc phase transition temperature con- 
centrational dependence was expressed through the transition temperatures of the composition and 
members. The obtained Ta {x) dependence was shown to fit pretty good the observed values for PZT 
and RADP. The application of the proposed theoretical approach to other types of mixed systems 
is discussed. 



I. INTRODUCTION 

Mixed antiferroelectric-ferroelectric systems like 
PbZri_,Ti,03 (PZT), Rbi_,(NH4),H2P04 (RADP) 
(0 < X < 1) attract much attention of scientists because 
of their unusual properties and phase diagram peculiar- 
ities. These anomalies made PZT material useful for 
many technical application [1]. In particular composi- 
tions of PZT enriched by Zr {x < 0, 1) appeared to be 
useful for application in the novel branches of electronic 
technique (see e.g. [2], [3] and ref. therein). Theoretical 
description of the properties and phase diagram nearby 
antiferroelectric phase transition is hard problem because 
of complexity of description of antiferroelectric phase it- 
self. The proximity of free energies of antiferroelectric 
and ferroelectric state is known to be the characteristic 
feature of antiferroelectrics [4]. This proximity lead to 
dielectric permittivity anomalies and to the possibility to 
induce ferroelectric order by application of small enough 
electric field. The latter can be the clue of the reasons 
of ferroelectric phase appearance even at small x values. 
Small content of ferroelectric material in the mixed sys- 
tem can be considered as the source of electric dipoles 
and induced by them random electric field. Calculation 
of this field magnitude allowing for dipole-dipole inter- 
action via ferroelectric and nonferroelectric soft modes 
had shown that mean field is induced by the interac- 
tion via ferroelectric soft mode whereas other soft modes 
increase mean square field [5]. The supposition that 
nonzero internal mean field can destroy antiferroelectric 
order by the same way as external field made it possible 



to describe semiquantitively the observed phase diagram 
[5]. 

In present work we performed detailed calculations of 
phase diagram of mixed antiferroelectric-ferroelectric sys- 
tem for the region of small concentration of ferroelectric 
material. The obtained data fitted pretty good measured 
phase diagram for PZT and RADP. The application of 
proposed method for calculation of phase diagrams in 
different mixed ferroelectrics is discussed. 



II. GENERAL EQUATIONS 

Let us begin with the description of antiferroelectric 
phase. In Ising model free energy F can be written as 
[4]: 
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In [2ch2P{p*E + Js<7) + 2ch2pjAS] 



(1) 



where /? = l/(fcT), effective dipole moment p* = (eq — 
l)jaP/3, Sa and ja are antiferroelectric component sus- 
ceptibility and Lorentz factor respectively, p = Zb (±6 
are the positions of two- well potential minima), E is 
electric field, a and S are respectively ferroelectric and 
antiferroelectric transition order parameters (ct = 
S = rj/b, ^ and rj being mean homogeneous and inho- 
mogeneous displacements respectively); J a = Ta and 
Js = Ts are respectively antiferroelectric transition and 
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Curie temperatures. Note, that in this model mean 
dipole moment of the unit cell has the form: 

dr^Z{l±ri) (2) 

Thus we consider the simplest two-sublattice model of 
antiferroelectric. 



A. Order parameters 



One can see that in mean field approximation when 
f{£,a,L) has the form of (5-function Eqs. (H) give Eqs. 
(||) at a; = and Eq. (|) at a; 1. 

Eqs. (|^) are general expressions for order parame- 
ters dependence on concentrations and characteristics of 
mixed system components. Note, that in the case S ~ Q 
(both components of mixed system are ferroelectrics) the 
equations for a and L have the same form as it has to be 
expected. 



Eq. (|l|) leads to the following expressions for the order 
parameters: 
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(3) 



ch2P{p*E + Jsa) + ch{2pjAS) 

Eqs. describe antiferroelectric component, i.e. 

mixed system at x = 0. Ferroelectric component corre- 
sponds to a; = 1 and its order parameter L in Ising model 
has the form [4] 



L = th(3{EoL + d*E) 



(4) 



where Eq — Tcmf is mean field, that defined ferroelectric 
transition temperature, d* is effective dipole moment of 
ferroelectric component ions. 

Addition of small amount of ferroelectric component to 
antiferroelectric one leads to appearance of random elec- 
tric field e, induced by electric dipoles of this component 
in the mixed system. Mean field in the mixed system has 
to include the contribution from ferroelectric component 
and from ferroelectric order parameter of antiferroelec- 
tric component namely 
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xEqL 
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This field as the most probable one has to define the 
position of maximum of random field distribution func- 
tion /(e, a, L). It is equal to 5{e — Js{l — x)a — xEqL) in 
mean field approximation. Beyond of this approximation 
one has to take into account random field, i.e. to substi- 
tute E + e {E ^ 0) OT e {E = 0) for E in Eqs. (|), (|) 
and to perform averaging with random field distribution 
function. This yields: 



sh{2/3p*{e + E))f{e,a,L)de 



ch(2(3p*{e + E)) + ch{2f3JA{l - x)S) 



B. Transition temperatures 

Let us proceed now to calculations of the transition 
temperatures Ta (to antiferroelectric phase) and Tc (to 
ferroelectric phase) dependence on concentration. For 
the phase transitions of the second order these temper- 
atures can be obtained from Eqs. (^) in the limit of 
L— >0. In such a limit one can simplify 
the integrands in Eqs. (||) substituting ch2(3 JaS « 1, 
sh2P JaS « 2/3 JaS and representing the distribution 
function in the form: 



/(£,a,i)=/o(e) 



./^(l — x)(7 EoxL\ / df 
de 



=L=0 

(7) 



Here /0(e) is the distribution function in the case of mean 
field absence. 

Keeping in mind that fo{^) is even parity function in 
aforementioned limit Eqs. (m give: 
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s 



sh{2pjA{l - x)S)f{e, a, L)de 
ch{2(3p*{e + E)) + ch{2f3JA{l - x)S) 



thij3d*{e + E))f{e,a,L)de 



Integral (H) has to depend on concentration because of 
dependence on x of the distribution function width. 
(6) The solution of the system of Eqs. (^) yields: 

r„ = (1 - x)JAla{x),T, = [(1 - x)Js + xEo] Ic{x) 

(10a) 

or in dimensionless variables 
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Eq. (10a) can be rewritten as 

Ta = (1 - x)Ia{x);Tc = [a;A + (1 - x)] Ic{x) 



(10b) 



(10c) 



Here Ia{x) and /c(a;) are integrals (g) at /3 = Pa and 
P = (3c respectively. One can see, that the transition 
temperatures are defined by the corresponding transi- 
tion temperatures and concentrations of the components 
of mixed system and random field distribution function 
characteristics. 

Equations (10) define concentrational dependence of 
transition temperature, i.e. phase diagram of mixed sys- 
tem for small concentration of ferroelectric component. 
One can see, that Eqs. (10) give the correct values in 
mean field approximation: 

Ta = Ta, Tc = Ts at X = and Ta = 0, Tc = i?o = 
Tcmf at a; = 1. 



C. Critical concentrations 

Random field induced by electric dipoles of ferro- 
electric component of mixed system can destroy com- 
pletely antiferroelectric phase at some concentration of 
the dipoles. This concentration x = Xca is known to be 
critical concentration at which Ta = 0, so that antiferro- 
electric phase exists at a; < Xca only. On the other hand 
ferroelectric phase in the mixed system arises at x > Xcc 
where Xcc is critical concentration for ferroelectric phase 
transition. It has to be Tc > Ta at a; > Xcc- Therefore 
critical concentrations Xca and Xcc in the region a; ^ 1 
can be obtained from the equations = {xca) and 
Tc ~ Ta {x — Xcc)- Let us limit ourselves by the case 
A = 0, which is of importance both for case of PZT and 
RADP on the " edges" of phase di agra m. To perform the 
calculations on the base of Eqs. (10c) we have to know 
concentrational dependence of integrals Ia,c{x) defined 
by the distribution function fo{y). For electric dipoles as 
random field sources it can be represented in Gaussian 
form (see [6] for details): 
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(12) 



Here eo, a and rc are respectively dielectric permit- 
tivity, lattice constant and correlation radius connected 
with soft ferroelectric mode of antiferroelectric compo- 
nent of mixed system. Parameter ^ > 1 takes into ac- 
count dipole-dipole interaction via antiferroelectric and 
other nonferroelectric soft modes [5]. Performing chang- 
ing of variables in integrals Ia,c {(3rnE = u), we obtain 



that this reduces just to renormalization of parameter cq 
in (^, giving Cm = com"^ , m = p*,d*. 

To map correctly the interval < z < oo into interval 
< a: < 1 we adopt following parametrization for Zc^a 



VcA^)P^^ P = —1 ^cix) = 
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Let us begin with critical concentration for antiferroelec- 
tric phase transiti on. 

At /3 = /3a Eqs.([lOa|) with respect to Eqs. (|ll|), ([l|) at 
Pa ^ oo yields 
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Sub stitu tion of Eq. (14) into the expression for (see 
Eqs. (lOc)) gives the following equation for critical con- 
centration Xca' 
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(15) 



Eq. ( p^ ) has been solved numerically during the calcula- 
tion of Ta- 

Critical concentration for ferroelectric phase transition 
can be obtained by the same way, but in the limit of 
Tc — > Ta, i.e. at Tc Ta/Ts, this ratio being larger or 
close to unity. Since in such a limit integral Ic{x) can be 
calculated only numerically we shall obtain Xcc value in 
the next section devoted to numerical calculations with 
two fitting parameters A = Tcmf/Ts and Qc = Cd*/Ts- 



III. PHASE DIAGRAM OF MIXED 
ANTIFERRO-FERROELECTRIC SYSTEM 

Phase diagram, i.e. concentrational dependence of an- 
tiferroelectric and ferroelectric phase transition temepra- 
ture was calculated with help of Eqs. (10). Integrals 
Ic.a{x) were calculated numerically. The obtained phase 
diagram is represented in Fig.l in the region of small con- 
centration of ferroelectric component for several values 
of dimensionless parameters Qa, Qc and A (see Eq.(lOb)). 



One can see from Fig.l that antiferroelectric phase transi- 
tion critical concentration Xca increases with qa decrease, 
its values being close to those calculated on the base of 
Eq. ( |l5|) . The obtained critical concentration of ferroelec- 
tric phase transition several times smaller than that for 
antiferroelectric transition. This seems to be the result 
of proximity of free energies and may be of Ts and Ta 
which is characteristic feature of antiferroelectrics. The 
dimensionless transition temperature to antiferroelectric 
phase Ta decreases with concentration increase for all the 
considered parameters, whereas the behaviour of that to 
ferroelectric phase Tc depends on the choise of the pa- 
rameters values. In particular Tc can increase or decrease 
with increase of electric dipoles concentration (see Fig.l). 
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Physically these types of behaviour depend on distribu- 
tion function width and ratio (A) of ferroelectric phase 
transition temepratures Tcmf of ferroelectric component 
of mixed system to that of antifcrroelectric component 
Ts- the smaller the distribution function width and the 
larger A value, the larger ferroelectric phase transition 
temperature in the mixed system. 



IV. DISCUSSION 

4.1. Let us begin with comparison of obtained data 
with experiment. The prominent example of mixed ferro- 
antiferro system is known to be PbZri_a;Tia;03 (0 < 
a; < 1), where PbZrOa is antifcrroelectric with transition 
temeprature from cubic paraelectric phase to rhombic an- 
tifcrroelectric phase at Ta = 500 K and PbTiOa is ferro- 
electric with transition temeprature from cubic paraelec- 
tric phase to tetragonal ferroelectric phase at Tcmf = 763 
K. Phase diagram of mixed PZT system was investigated 
in many papers, but mainly at T > Troom (see e.g. [6] 
and ref. therein). Therefore critical concentration Xca 
is known, and = 300 K at x = (5 — 6)%. One can 
see that the curve for = 1.7 gives Ta ~ 300 K at 
X = 0.06 and Xca ~ 0.11. General form of calculated Tq 
concentrational dependence looks like experimental curve 
at X < 0.1 ( Fig.l). There is a scattering of experimen- 
tal data for ferroelectric phase transition temperature at 
small PbTiOa content {x < 0.05) because of influence 
of some uncontrolled defects and impurities on Tc value. 
The most probable behaviour of Tc is its slight increase 
(Fig.l) [6]. Such behaviour is similar to the calculated 
curve with qc = 1.7 and A = 7. 

Another important example of mixed ferro-antiferro 
system is Rbi_:r(NH4)a;H2P04 (0 < a; < 1) as a repre- 
sentative of hydrogen-bonded material. It was studied 
in many details (see e.g. [7]). RbH2P04 (RDP) under- 
goes a ferroelectric phase transition at Tcmf = 146 K, 
NH4H2PO4 (ADP) becomes antifcrroelectric near Ta = 
148 K. In mixed RADP system antifcrroelectric transi- 
tion temperature Ta decreases with increase of RDP con- 
tent and Tq = K near x = 0.2 of RDP admixture, i.e. 
Xca ~ 0.2 [8]. Substitution of this Xca value into Eq.(p^) 
gives Qa ~ I. One can see from Fig. 2 that the curve, 
which corresponds to qa = 1.175 fits the experimental 
points. Contrary to PZT there is no data about exis- 
tence of ferroelectric phase transition at small (< 10%) 
content of RDP. The reason for this difference can be the 
approximately equal values of Tcmf and Ta in RADP 
whereas they are different in PZT. Because of this the- 
oretical parameter A « Ta/Ts ~ 1 in RADP and at 
A = 1 Tc = Tc/Ts curve may be lower than Ta curve and 
only antifcrroelectric phase transition exists. One can see 
from Fig. 2, that proposed theory explaines the observed 
behaviour of ferroelectric phase transition decrease with 
increase of antifcrroelectric component concentration and 
observed critical concentration Xcc {qc = 1.424, A = 1.1). 



Note that tuneling has to be taken into account for mixed 
systems of KDP family. Preliminary calculations have 
shown, that tuneling result into change of fitting param- 
eters values. 

Therefore the proposed theory describes pretty good 
phase diagrams of PZT and RADP for small x values. 
Note, that only for small x {x < 0.2) Ti and Rb ions 
can be considered as impurities randomly distributed in 
PbZrOs and NH4H2PO4 matrices so that they can be 
the sources of random field. 

4.2. The proposed theory can be applied for another 
types of mixed systems. For example, mixture of two 
ferroelectric systems {Fi)i-xiF2)x at a; < 0, 1 described 
with the help of Eqs. (|), (|) at zero antifcrroelectric 
phase order parameter (S — 0). In this case the Equa- 
tions for a and L completely the same, i.e. a = L. 
In such situation mean field = Tcmf is ferroelectric 
phase transition temperature of Fi component whereas 
F2 component defines the form and width of random 
field distribution function. Because the main physical 
idea of the proposed model is the inhibition and distruc- 
tion of long range ferroelectric or antifcrroelectric order 
by random electric field of different sources this model 
can be useful for description to relaxor ferroelectrics 
like PbMgi/3Nb2/303 (PMN), PbSci/2Tai/203 (PST), 
Pbi_yLaj,Zro^65Tio,35 3 (PLZT, y = 0,08; 0,09). But 
the relaxors like PMN, PST are not solid solutions of 
two component, and in PLZT La ions arc the impurities. 
Thus in all these materials there is onlu one ferroelectric 
component known as Burns reference phase (see [9] and 
ref. therein), which is PbZro,65Tio.35 3 in the of PLZT. 
The transition temperature from para- to ferroelectric 
phase of the reference phase Tcmf = defines the po- 
sition of random field distribution maxima whereas the 
substitutional disorder, vacancies of lead and oxygen, im- 
purities like La in PLZT define the distribution function 
width. Namely such a model has been proposed recently 
[10]. It permitted to describe the observed phase dia- 
gram of the relaxors, their dynamic and static properties 
peculiarities (see e.g. [11] and ref. therein). 

Since the model is able to describe the distraction of 
long range order by any random field sources it can be 
useful also for description of solid solutions of ferroelec- 
tric and dielectric components, e.g. Ba(Tii_2:Sn2;)03 
{x < 0,1). The observed in this material relaxor-like 
behaviour as well as change of phase transition order [12] 
can be the consequence of strong random field infiuence 
as it was shown recently [11,13]. Note that some residual 
ferroelectric domains were observed in another relaxor 
PLZT (9/65/35) [14]. As the matter of fact, the dipole 
glass and mixed ferroglass phases (which appear in the 
relaxors with temperature lowering) can be considered as 
the reentrant phases with respect to this residual ferro- 
electric phase existing near the Burns temperature. 

Therefore the generality of proposed model gives a pos- 
sibility to apply it to various mixed systems which are 
solid solutions of ferro- and antiferroelectrics, two ferro- 
electrics, ferroelectric and dielectric etc. as well as to 
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ferroelectric relaxors. 
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